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INTRODUCTION 


Modem  tank  cannons  are  long,  relatively,  thin,  beam-like  hollow  cylinders.  Their 
accuracy  is,  in  part,  determined  by  their  flexibility,  especially  under  dynamic  loading.  Very 
small  deflections  and  rotations  of  the  muzzle  end  can  have  a  significant  influence  on  the 
accuracy  of  the  shot  at  long  ranges.  Muzzle  motions  induced  by  firing  are  inevitable,  and 
difficult  to  control  because  the  time  scale  of  firing  is  of  the  order  of  milliseconds. 

Another  source  of  muzzle  motion  is  the  ground-induced  motion  of  the  vehicle.  These 
motions,  transmitted  through  the  trunnions  and  gun  actuators,  can  be  quite  large  and  have 
frequencies  comparable  to  the  natural  frequencies  of  the  tube.  The  time  scales  of  these 
disturbances  depend  on  the  tank  speed  and  on  the  nature  of  the  terrain.  They  are  typically  of  the 
order  of  seconds  or  longer.  Sensing  and  actuation  to  control  the  influence  of  vehicle  motion  on 
the  muzzle  response  might  be  possible,  given  these  relatively  long  time  scales.  This  raises  two 
questions.  First,  is  it  possible  to  reject  some,  or  all  ground  motion  disturbance  from  the  muzzle 
motion?  In  a  previous  paper  (ref  1)  it  was  suggested  that  not  all  of  the  disturbance  could  be 
rejected.  Second,  if  the  more  comprehensive  model  used  here  indicates  that  all  of  the 
disturbance  can  be  rejected,  what  is  the  required  control  strategy? 

During  a  discussion  with  Dr.  Purdy,  author  of  Reference  2,  he  suggested  that  the  fidelity 
of  his  two-segment  flexible  model,  documented  in  Reference  1,  was  inadequate.  He 
recommended  that  the  tube  should  be  divided  into  at  least  three  segments,  with  intervening 
torsional  springs  and  dampers.  The  author  is  indebted  to  Dr.  Purdy  for  this  suggestion  since  this 
report  is  the  result  of  his  recommendation. 

EQUATIONS  OF  MOTION 

Figure  1  shows  the  generic  model  of  the  tube  and  the  various  quantities  that  determine  its 
dynamic  behavior. 


Figure  1 .  Gun  generic  model. 
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A  free  body  analysis  of  this  model  yields  the  classical  dynamic  equation 


HPl+lW+lM-M  (i) 

The  elements  of  mass,  damping,  and  stiffness  3x3  matrices  are  shown  in  Appendix  A. 

Purdy  (ref  3)  has  shown  that  the  tube  motion  can  be  adequately  modeled  if  the  segmented 
model  matches  the  pinned-free  and  cantilever  frequencies  of  the  mounted  tube.  Matching  is 
accomplished  by  adjusting  the  size  of  the  rigid  segments  and  the  stiffness  of  their  connecting 
torsional  springs.  The  2x2  matrices  for  the  cantilever  mode  are  outlined  in  Appendix  B. 

Transformation  of  equation  (1)  into  the  frequency  domain  will  allow  its  incorporation 
into  the  control  strategy.  Taking  the  Laplace  transform  of  equation  (1)  yields 


(2) 


The  elements  of  [a]  and  [7]  are  listed  in  Appendix  C.  This  equation  relates  the  response 
vector  on  the  left  to  the  disturbance  vector  on  the  right.  These  vectors  also  contain  the  actuator 
force,  Fc,  and  the  actuator  displacement,  xc,  in  addition  to  disturbances  (s2y,(s),  s8p(s),  0p(s))  and 
the  responses  0j(s),  Qi(s),  Qi(s). 


*,(*) 

F.(f)  ' 

0,(j) 

-=HWW' 

s‘y,(s) 

s6,(s) 

M), 

*,«  J 

(3) 


where  [C]  is  the  transpose  of  the  numerators  of  the  cofactors  of  [a].  The  elements  of  [C\  and  det 
[a]  are  listed  in  Appendix  D. 

The  final  step  in  the  preparation  of  the  dynamic  equations  is  to  perform  the  operation 

lk\[cWMB]  (4) 

where  the  elements  of  [5]  can  be  found  in  Appendix  E. 
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The  result  of  these  straightforward,  but  laborious  manipulations,  is  an  equation  for  the 
response  to  the  disturbance  in  terms  of  the  properties  of  the  model  contained  in  [5],  i.e., 


(5) 


Of  course,  it  was  known  at  the  outset  that  equation  (1)  could  be  put  into  this  form.  This 
section  merely  provides  the  details  of  how  this  transformation  is  performed,  and  documents  the 
intermediate  steps  and  their  components. 

FEEDBACK  AND  FEEDFORWARD  CONTROL 


The  portion  of  the  response  due  to  the  applied  actuating  force  is 


Xc(s) 

'V 

0,(s) 

*21 

02(s) 

*3, 

Ms). 

*41. 

(6) 


where  [Gp]  is  the  "plant"  transfer  function. 

The  portion  of  the  response  due  to  the  disturbance  is 


xM 

'*.2 

Bn 

B„ 

0,(s) 

b2Z 

Bn 

*24 

B32 

Bn 

*34 

p 

B,2 

Bn 

*44. 

s2y,(s) 

sOp(s) 


(7) 


where  [Gd]  is  the  disturbance  transfer  function  and  {D}  is  the  disturbance  vector. 

Figure  2  is  a  block  diagram  of  the  gun  tube  system  with  a  gain  Gc,  feedback  H, 
feedforward  GCd,  and  reference  signal  R.  Because  R  is  a  scalar,  the  feedforward  transfer  function 
is  a  row  vector,  i.e., 


[gJ=[g„g„g„] 


(8) 


Figure  2.  Block  diagram  of  gun  tube  system. 


Assuming  that  the  tube  rotations  at  the  trunnions  and  the  muzzle  can  be  sensed,  the 
feedback  transfer  function  is  also  a  row  vector,  i.e., 

M  =  0  G„]  (9) 


Now  referring  to  Figure  2,  the  response  to  the  disturbance  for  the  controlled  system  is 


N-lG,J[G,][G0,lM  =  I/]+[G,J[<?e][/f]|fe,}  (10a) 


or  in  abbreviated  notation 


(10b) 


The  final  step  in  these  manipulations  is  to  solve  for  the  response  to  the  disturbance,  which 
is 


K } = W'  M  {£>} = ^  W  M  M  doc) 

where  [r]  is  the  matrix  of  the  co factors  of  [<y]. 

The  matrices  [r]T  and  [d\  are  given  in  Appendix  F.  It  is  interesting  to  note  that  [d] 
contains  only  the  feedforward  transfer  functions  Gn,  Gn ,  and  Gn,  while  [r]  and  de([q]  contain 
only  the  feedback  transfer  functions,  G22,  and  G24. 
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DISTURBANCE  REJECTION 

Referring  to  Appendix  F,  the  expanded  version  of  equation  (10c)  is 

\ril  r21  0  r4, 

0/M  1  0  r22  0  r42 

det[q\  0  r23  r33  r4i 

?i(s)\D  1°  r24  0  r44 

where 

det[q\  =  1  +  B2i  Gc  G22  +  B41  Gc  G24 

To  completely  remove  the  effect  of  the  disturbances  on  Os(s)d  requires  that 

r24  d 21  +  r44  d4i  ~  0 
r24  d22  +  r44  d42  ~  0 
r24  d  23  +  r44  d43  =  0 

and 

det\q\  -*■  0 

One  way  to  accomplish  this  is  to  let  G22  =  0  so  that 

r24  =  ~  B41  Gc  G22  =  0 

and  then  choose 

d4l  =  B42  —  B4]  Gc  Gn  —  0 
d42  =  B43  —  B4I  Gc  G,2  =  0 
d43  =  B44  —  B41  Gc  Gi3  =  0 

so  that 

det[q\  =  1  +  B4]  Gc  G24  =  rn  =  r22  =  r33 


(lOd) 

(11) 

(12a) 

(12b) 

(12c) 

(12d) 

(13) 

(14a) 

(14b) 

(14c) 

(14d) 
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and 


r44  ~  1  ’  r21  ~  r23  ~  r24  ~  ^ 

The  effect  of  this  choice  on  the  disturbance  transfer  function  is 


1 

det[q ] 


WM-— 


1 

0 

0 

£ 
_ I 

'dn 

dl2 

d13' 

1 

0 

r22 

0 

r42 

d2, 

d22 

d23 

r„ 

0 

0 

r33 

r43 

d3I 

d32 

d33 

0 

0 

0 

r<4_ 

_  0 

0 

0 

(15) 


(16) 


An  alternative  strategy  is  to  let  r44  =  1  +  Bji  Gc  G22  =  0,  and  then  choose  d 21  =  d.22  =  <4?  = 
0  so  that  det[q]  =  B4j  Gc  G24  =  rn  -  r33.  However,  in  the  end  the  resulting  disturbance  transfer 
function  is  the  same  as  created  by  equation  (16). 

Since  the  fourth  column  of  [r]ris  eliminated  by  the  matrix  multiplication  and  det[q\  =  rit 
=  *22  =  r33,  det[q]  will  be  eliminated  from  the  transformation. 

IMPLEMENTATION 

This  analysis  of  the  three-segment  model  indicates  that  model  muzzle  element  can  be 
stabilized  by  properly  selecting  the  feedforward  and  feedback  transfer  functions.  This  is  contrary 
to  the  finding  for  the  lower  order  two-segment  model  (ref  1).  Although  the  three-segment  model 
only  approximates  the  real  tube,  the  results  of  this  model  are  encouraging  with  respect  to  real 
tubes. 


In  order  to  achieve  muzzle  stabilization,  the  breech-end  of  the  gun  must  be  actuated. 
Segments  1  and  2  will  also  rotate.  These  motions  are  determined  by  equations  (14)  and  (16).  All 
of  the  elements  of  [d\  can  be  written  in  terms  of  the  elements  of  [5].  Because  [r]r  acts  like  an 
identity  matrix,  the  product  [r]r[<f]  is  quite  simple.  If  the  process  is  carried  a  step  further,  the 
result  can  be  put  in  terms  of  the  elements  of  [C]  with  startling  results,  i.e.,  *4;  =  c4?  =  d32  =  d33  = 
0,  and  only  dn  -  dn  -  d/3  =  d2i,  and  d3i  are  nonzero.  The  surviving  elements  are 

d n  —  —  {{^14  C31  ~  C/i  C34 )  1 32  +  {C/4  C41  ~  C/i  C44  )l42 )  (17a) 

-  c  -d-trs  ((C„  c„  -  C„  cu)lu+  (C„  C4}  -  C„  Cjlj  (17b) 

dj  1  =  ((^/.  CJ}  -C„C„)IU  +  (CN  C43  —  C„  Cjj  )  /., )  (17c) 
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(17d) 


d „  -  -  j  ,  1  (c„  c„  -  C„  Cj/„  (17e) 

These  transfer  functions  relate  the  disturbances  to  the  responses.  All  of  the  Cs  are  of 
order  s4  with  the  exception  of  Cu,  which  is  of  order  s6.  Since  there  is  no  restraining  torsional 
spring  connecting  the  tube  to  the  mount  in  the  model,  s  =  0  is  a  root  of  det[a\.  Removing  this 
rigid  body  factor  from  det[a ]  reduces  it  to  order  s5. 

TUBE  MODEL  PARAMETERS 

The  feedforward  transfer  functions  depend  on  the  length  and  mass  properties  of  the 
segments,  the  torsional  stiffness  of  the  joining  springs,  and  the  torsional  damping  coefficients. 
These  are  chosen  so  that  the  actual  cantilever  and  pin-free  mode  shapes  and  natural  frequencies 
are  matched  as  closely  as  possible  (ref  3).  To  simplify  the  matching  process,  it  is  assumed  that 
the  damping  is  negligible.  The  first  estimate  of  the  segment  lengths  can  be  obtained  by  "fitting" 
the  straight-line  segments  to  the  mode  shapes  obtained  from  a  finite  element  model  of  the  tube  or 
other  modal  analyses.  This  fitting  is  best  done  by  graphically  overlaying  the  segments  on  plots 
of  the  mode  shapes  to  estimate  the  segment  lengths.  The  calculation  of  the  mass  properties  of 
the  segments  can  then  be  performed  and  these,  along  with  the  modal  frequencies,  inserted  into 
the  characteristic  equations.  The  characteristic  equations  will  then  contain  only  the  torsional 
stiffnesses  as  unknowns.  The  cantilever  and  pin-free  equations  are  both  quadratic  so  that  the 
stiffness  coefficients  can  be  found  directly.  The  degree  of  matching  is  determined  by  how 
closely  the  cantilever  and  pin-free  stiffnesses  agree. 

The  characteristic  equations  for  the  cantilever  and  pin-free  segments  are  given  in 
Appendix  G.  Although  the  pin-free  equation  appears  to  be  sixth-order,  it  has  a  double  root  that 
is  zero.  The  calculations  for  this  trial-and-success  process  are  easily  implemented  on  a 
spreadsheet. 

The  XM291  tank  gun  was  chosen  for  modeling  because  its  mode  shapes  and  frequencies 
were  available  from  an  existing,  validated  analytical  model.  Matching  the  stiffnesses  proved  to 
be  surprisingly  easy,  requiring  only  modest  adjustments  to  the  first  estimates  of  the  segment 
lengths.  Since  all  their  frequencies  (cantilever:  97.4  Hz,  40.35  Hz;  pin-free:  25.08  Hz,  81.59  Hz) 
were  inserted  into  the  characteristic  equations,  they  are  matched  exactly.  The  torsional 
stiffnesses  for  the  pinned- free  and  cantilever  modes  were  matched  within  2%  using  the  lengths  // 
=  6.0  ft,  h  =  5.5  ft,  and  h  -  6.0  ft.  From  this  process  the  model  torsional  stiffnesses,  kn  = 
3.6(106)  lb-ft/rad  and  k23  =  1.69(106)  lb-ft/rad  were  obtained 

Dynamic  analyses  (refs  2,  3)  have  successfully  modeled  tube  response  using  proportional 
damping,  i.e.,  [c]  =  f^k].  In  the  case  of  the  XM291,  /?=  0.0015  second  has  been  found  to  be 
reasonable.  A  reasonable  estimate  for  trunnion  damping  is  cIp=  750  lb-ft-s/rad. 
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The  first  attempt  to  determine  the  feedforward  transfer  function  using  equations  (12)  and 
(13)  failed  because  some  of  the  roots  of  B4i  were  positive.  This  difficulty  was  eliminated  by 
using  the  alternative  strategy  described  above,  with  the  following  results: 


r44  =  1 +  B2,  Gc  G22  =  0  (18) 

G.  G„  =  —  (19a) 

G,G„=^  09b) 

Gc  G,3  =  (19c) 

B2I 


The  det[a ]  plays  no  role  in  these  fimctions  because  it  is  canceled  by  ratioing  the  B's. 
Figures  3  and  4  show  the  Bode  plots  of  Gc  Gu  and  Gc  G12 ■  The  transfer  function  Gc  G13  is  zero 
so  that  0P  is  not  fed  forward.  The  numerators  and  denominators  are  all  fifth-order  polynomials, 
so  that  the  high  and  low  frequency  gains  are  bounded. 


Frequency  (rad/sec) 


Figure  3.  Bode  plots  of  Gc  Gu. 


8 


Phase  (deg):  Magnitude  (dB) 


Figure  5.  Bode  plots  of  Gc  G22. 
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The  elements  of  d22,  d2s,  d22,  and  da  of  [d\  were  found  to  be  identically  zero.  The 
remaining  nonzero  elements  of  equations  (lOd)  and  (16)  yield  the  following  response  equations: 


>,{s)+dl2sOp{s)  +  d13  9p(s) 

(20a) 

0,  =d2,s2  y,{s) 

(20b) 

02  =d3Is2y,(s) 

(20c) 

Figures  6  through  10  show  the  transfer  functions  required  by  the  equations  above. 
Figures  9  and  10  show  that  the  effects  of  the  trunnion  acceleration  on  &i  and  02  are  highly 
attenuated  so  that  large  angular  displacements  of  the  tube  are  not  required  to  achieve 
stabilization. 

Figures  6  through  8  are  quite  similar.  It  appears  that  the  required  xc  will  depend  largely 
on  the  trunnion  acceleration  and  pitch  rate  at  very  low  frequencies.  There  is  a  considerable 
attenuation  of  the  disturbance  inputs  up  to  100  rad/second  (~15  Hz)  with  a  return  to  the  low 
frequency  levels  at  103  rad/second  (160  Hz). 


Frequency  (rad/sec) 


Figure  6.  Bode  plots  of  du. 


Phase  (deg);  Magnitude  (dB)  Phase  (deg);  Magnitude  (dB) 
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Frequency  (rad/sec) 

Figure  7.  Bode  plots  of  dj2- 


Frequency  (rad/sec) 

Figure  8.  Bode  plots  of  du¬ 


ll 


Phase  (dag);  Magnitude  (dB)  Phase  (deg);  Magnitude  (dB) 


CONCLUSIONS 


It  appears  that  the  results  previously  obtained  with  the  two-segment  model,  (ref  1),  led  to 
the  erroneous  conclusion  that  the  effects  of  the  disturbance  could  not  be  entirely  repeated  from 
the  muzzle  angular  displacement.  The  analysis  of  the  three-segment  model  presented  suggests 
that  this  is  possible,  at  least  theoretically.  Of  course,  the  unanswered  question  is,  "What  would 
be  revealed  by  a  higher-order  multi-segmented  model,  and  how  many  segments  are  enough?" 

On  the  practical  side,  it  is  certain  that  the  transfer  functions  cannot  be  duplicated 
precisely.  There  are  four  of  these  that  must  be  implemented  with  reasonable  fidelity  to  achieve 
the  predicted  results  of  the  three-segment  model.  That  number,  along  with  the  input  signals, 
indicates  the  magnitude  of  this  task.  While  feedforward  and  feedback  control  have  long  been 
used  in  fire  control,  it  is  hoped  that  this  report  provides  some  guidance  in  their  use  when  tube 
flexure  is  a  consideration. 
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APPENDIX  A 


[m]  = 


mi2 

m,3 

m2, 

m22 

m23 

m3, 

m32 

m33 

mu  =  /,  +  m{r]\  +  lf(m2  +  m 

m22  ~  1 2  ^2^2  ^3^2 

w33  =  /3  +  mjil 

mi2  =  m2\  =  m2^\rh  +  W3/,/2 

wn  =  m31  =  m3/j?73 

w23  =  myi  =  w3/2^3 


M 


—  cn  0 

C12  +  C23  ~  C23 

—  c23  c23 


k  12 

kj2  +  k23 


-k 
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APPENDIX  B 


m22 

m32 


m33 


C12  +  C23 


'23 


kj2  +  k23 

-k23 

—  k23 

k  23 

APPENDIX  C 


H- 


a,,  al2 

a2i  a22 


0  a 
0  a 


32 


42 


0  0 


a23  a24 

a33  a34 

a  43  a44 


a„  =  kd 

an  =  a  j,  =  -KXp 


a24  ~  a42  ~  mi3S 


CI22  —  MjjS  +  (cI2  +  cIp  )s  +  (kdXtp  +k12) 
a  23  =  o.32  =  mI2s  —  cl2s  —  kl2 
ci  33  —  th22s  +  ( c!2  +  c23  )s  +  ( kI2  +  k23 ) 
a  33  ci 33  —  m23s  c 2s s  k  23 

Q44  ~  M33S  ^23^  k23 
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0 


0 


[/]= 


0 

0 

0 


1 22 


32 


42 


0  0 
0  0 


In  =  1-0 
1,4—KX* 


'22 


1  23 


1  24 


1  32 


42 


=  ~( V,m,  +l,m2  +  l,m3) 

=  c>p 

=  ~(V2m2+12m3) 

=  ~rl3m3 
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APPENDIX  D 


o' 

C,2 

C,3 

C,4 

C, 

C22 

C23 

c24 

'Cu 

c32 

C33 

C34 

C4, 

C<2 

c43 

c44 

('ll  =  a22a33a44  +  2a23a34°24  ~(a24a33°42  +  a43a34°22  +  a23°32a44  ) 
(-22  ~  ail(a33a44  ~a34a43 ) 

(■33  =  ail(a22a44  ~a24a42) 

( 44  ~  atl(a22a33  ~a23a32  )~ai2a2!a33 
('12  =  ('21  =~°2l(a33a44  ~  a34°43  ) 

C 13  —  CSI  =  Q2l(a33a44  ~a34a42  ) 

C i4  —  C4I  =  —a2l(  a32a43  —  a33a42 ) 

( 23  =  C32  ~  ~au( a32a44  ~a34a42  ) 


('24  ~  ('42  ~  ail(a32a43  a33°42  ) 

( 34  ~('43  =  ~ail( a22a43  ~  ° 23° 42  )  +  °  12° 21°43 

det\a]  =  a,i(a22a33a44  +  a23a34a42  +  a32a43a24 ) 

~  ail(a24a33a42  +  a34a43a22  +  0 23° 32° 44  )  ~  a  12° 2l( 0 33a 44  ~  a34°43  ) 
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APPENDIX  E 


M- 


Bn  BI2  Bi3  BI4 
B 21  B22  B23  B24 
B31  B32  B33  B34 
B41  B42  B43  B44 


Bn  —  C„In  /  det[a\ 

B 12  =  (C2J22  +  C3II32  +  C4II 42 )  /  det\a\ 
B  a  =(C2lI23)/det[a\ 

B 14  =(C,,II4  +C2II24)/det[a\ 

B  2j  +(Cl2Iu)/det[a] 

B22  ~  (C2J22  +  ^32^32  +  C 42 1 42  )  /  det[o] 
B  23  =  (C22I23)/det[a\ 

B24  =  (CI2IJ4  +  C22I24 )/  det\a\ 

B3,  =(C,3I„)/det[a] 

B32  ~  (CJ22  +  C 33 1 32  +  C43I42  )  /  det[a\ 
B33  =(C23I23)/det[a] 

B34  —  (C,3I,4  +  C 23 1 24  )  /  det\a\ 

B  41  =(C,4I„)/det[a\ 

B  42  —  (C 24 1 22  +  C 34 1 32  +  C44I42 )  /  det[a\ 
B43  =  (C24I23)/det[a] 

B44  =(CI4Il4  +C24I24)/det[a] 


APPENDIX  F 
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0 

1 - 

0 

r22 

0 

'42 

0 

r23 

r33 

'43 

0 

r24 

0 

'*44  _ 
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ri  i  —  det\a\  —  1  +  B2tGcGn  +  B4IGCG24  r33  —  ru 


r2i  =-b„gcg22 

r22  =  1  +  B4,GcG24 
r23  ~  ~B3IGCG22 
r24  =-B4,GcG22 


d„ 

d12 

d13 

M- 

d  21 

d22 

d  23 

d3, 

d  32 

d33 

d41 

d  42 

d  43  _ 

d,2  = 

Bis- 

~Bn 

GcG12 

d22  = 

B23- 

-B:, 

GcG,2 

d}2  -  B33  B3,GCGI2 
d42  =  B4 3  -  B4IGcG,2 


r4t  =  -B,,GcG24 
r42  =  ~B  2  ,GCG  24 
r43  =  -B3/GcG24 
r44  =  1  +  B2!GcG22 


du  =BI2  -B„GcGn 
d2I  =  B22  -  B2!GcGii 
d3i  =  B32  -B3iGcGu 
d 4i  =B42  -B4,GcGn 
d  13  =B14  -BnGcG,3 
d 23  ~B24  —  B2iGcGi3 
d3  3  =B34  B  3jGcG  j3 

d 43  =B44  -B4iGcGi3 


APPENDIX  G 

m,  =tu2  =  first  two  natural  frequencies 

kn  ~  h7!  +  &22  \m22m33  -  m23m32 )/(m3S  +  2m2i  +  m22 

+  mJ3 (g 7/ \m22m33  - m23m32 )/(ms3  +  2m 23  +m22)=0 

k,2  =  +  &2  \m22m3J  -  m23m32  )/m33-  (m33  +  2m2i  +  m22  ]k23  /  m33 
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NOMENCLATURE 


[a] 

[B] 

[c] 

Cl2,  C23 
C]p 

[Q 

[d\ 

{D} 

Fc 

Gc 

GCd 

Gd 

GP 

Gji,  G12,  Gu 
G22,  G24 
H 

M 

m 

kl2,  k23 
kd 

h,  h ,  h 

[w] 


Dynamic  matrix 

IQU\ 

Damping  matrix 
Damping  coefficients 
Trunnion  viscous  friction  coefficient 
Cofactor  matrix  of  [a] 

Disturbance  input  matrix 
Disturbance  vector 
Elevation  actuating  force 
Scalar  gain 

Feedforward  transfer  function  vector 
Disturbance  transfer  function 
Plant  transfer  function 

Feedforward  transfer  function  vector  components 

Feedback  transfer  function  vector  components 

Feedback  transfer  function  vector 

Forcing  function  matrix 

Stiffness  matrix 

Stiffness  coefficients 

Drive  line  stiffness 

Segment  lengths 

Mass  matrix 
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Disturbance  response  matrix 
Elevation  actuator  displacement 
Distance  from  trunnions  to  drive 
Vertical  displacement  of  the  trunnion 
Cofactor  matrix  of  [q] 

Proportional  damping  coefficient 


Center  of  mass  coordinates 
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